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1 . Introduction  and  Summary 

The  Investigation  reported  here  was  undertaken  to  determine  If  a 
stochastic  Interrelationship  exists  between  the  running  times  and  the  down 
times  of  the  Robinson  Nuclear  Power  Plant.  The  data  was  provided  to  us  by 
the  Probabilistic  Analysis  Staff  of  the  Nuclear  Regulatory  Commission  (NRG). 
Our  goal  was  to  understand  what  the  data  were  telling  us  about  the  rela- 
tionship between  the  series  of  running  times  and  the  series  of  down  times. 
One  way  of  achieving  this  goal  Is  to  obtain  a Box-Jenklns  (1976)  type  of 
"transfer  function  model"  between  the  running  times  and  the  down  times. 

The  Interpretation  and  uses  of  a transfer  function  model  for  the  situation 
considered  here  are  discussed  In  Section  1.1. 

A first  step  In  the  analysis  of  the  data  was  Its  careful  screening. 
This  was  done  In  order  to  eliminate,  those  observations  that  were  judged  to 
be  questionable  or  that  had  arisen  under  unusual  circumstances.  Such  ob- 
servations Introduce  spurious  autocorrelations  and  cross  correlations, 
and  thus  tend  to  obscure  the  Identification  of  a simple  relationship  that 
may  exist  between  the  running  times  and  the  down  times. 

We  would  like  to  emphasize  that,  for  an  analysis  of  data  of  the  type 
discussed  here  (often  referred  to  as  "messy  data"),  an  examination  and 
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screening  of  the  data  prior  to  model  building  are  Important  preliminary 
operations.  If  one  neglects  to  perform  these  operations,  one  may  face 
the  frustrating  task  of  attempting  to  fit  several  transfer  function  models, 
none  of  which  may  be  satisfactory. 

In  Figures  1.1  and  1.2,  we  display  a time  sequence  plot  of  the 
screened  down  times  and  the  corresponding  screened  running  times 

, t=l,2,...,28  . Note  that  represents  the  first  down  time,  X2 

the  second  down  time,  and  so  on,  whereas  denotes  the  first  running 

time,  Y^  the  second  running  time,  and  so  on.  Note  also  that  the  two 

plots  are  not  drawn  to  the  same  scale.  In  Figure  1.3  we  indicate  the  rel- 
ative positions  of  the  the  > t=l,2,...  . 
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Figure  1.3  State  of  the  system  versus  time. 


In  transfer  function  model  building  observations  must  be  considered 
In  pairs.  In  our  case,  the  paired  observations  will  be 

(X^,Y^),  (X2,Y2),  ...,  (X28.Y28)  . 

In  considering  the  above  pairs,  we  will  have  to  bear  in  mind  that  X^  pre- 
cedes Y^  In  chronological  time. 

In  Table  1.1  we  present  the  actual  values  of  the  screened  down  times 
X^  and  the  corresponding  screened  running  times  Y^  , t=l,2,...,28  . 

If  changes  in  a series  of  observations  Y^  , t=l,2,...  , tend  to  be 
anticipated  by  changes  in  another  series  of  observations,  say  X^  , t*=l,2,... 

- 2 - 
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then  is  said  to  be  a leading  indicator  of  . In  our  case  it  is 

reasonable  to  assume  that  the  down  times  are  a leading  indicator  of 

the  running  times  Y^  . Based  upon  this  consideration,  and  together  with 

an  analysis  of  the  available  data,  albeit  an  insufficient  amount,  the  best 
transfer  function  model  we  have  identified  and  fitted  is  given  as 

Y^  - 10.630  = 6.6X  - 0.55X  , . (1.1) 

t t t-1 


1.1  Interpretation  and  uses  of  the  transfer  function  model 

Transfer  function  models  are  generally  used  to  forecast  the  future 
values  of  a time  series  Y^  (in  our  case  the  running  times)  given  the  pre- 
vious values  of  the  leading  Indicator  series  X^  (the  down  times  in  our 

case).  However,  forecasts  of  the  running  times  based  upon  the  previous  and 
present  values  of  the  down  times,  via  Equation  (1.1),  must  be  used  with 
great  caution  for  the  following  two  reasons.  First,  the  occurrence  of  un- 
foreseen but  rare  circumstances  may  cause  the  future  running  times  to  be 
abnormally  small  (or  even  large).  An  example  of  this  is  a reactor  shut- 
down due  to  an  unforeseen  operator  error.  Second,  since  Equation  (1.1) 
is  based  on  only  28  observations,  it  may  not  be  too  reliable  as  a model 
for  forecasting.  We  can,  however,  make  several  observations  of  practical 
interest  based  upon  Equation  (1.1). 

We  first  note  that  previous  values  of  Y^  , such  as  2’ ' ‘ ’ ’ 

etc.,  do  not  appear  in  Equation  (1.1).  This  implies  that  the  running  time 
history  gives  us  little  information  about  the  individual  future  running 
times;  that  is,  the  next  running  time  may  be  unpredictable  from  a knowledge 
of  the  previous  running  times.  However,  future  running  times  on  the  average 
may  be  estimated  from  previous  running  times. 

An  Important  consequence  of  Equation  (1.1)  is  that  the  running  times 
Y^  appear  to  be  strongly  influenced  by  the  immediately  preceding  down 

times  X^  . Since  the  down  times  generally  correspond  to  maintenance  actions, 

we  can  make  the  following  conjecture: 

- 6 - 
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Conjecture:  Barring  unforeseeable  circumstances,  and  confining  ourselves 

to  the  limits  of  the  observed  data,  the  operating  times  are,  on  the  aver- 
age, increased  by  a factor  of  about  six  per  unit  increase  in  the  mainte- 
nance times. 

An  explanation  to  support  the  above  conjecture  is  that  the  more 
thorough  a job  of  repair  that  is  performed,  the  longer  the  next  running 
time  becomes.  This  is  perhaps  one  of  the  most  important  conclusions  that 
can  be  reached  from  our  analyses. 

Since  the  coefficient  of  -0.55  associated  with  X , is  small 

t-1 

compared  to  the  coefficient  of  6.6  associated  with  , we  will  ignore 

the  effect  of  • Even  though  the  transfer  function  model  is 

obtained  after  an  Involved  analysis,  as  discussed  in  the  remainder  of  this 
report,  the  simplicity  of  Equation  (1.1)  suggests  that  a plot  of  versus 

X^  , t=l,2,...,28  , should  be  approximately  linear.  The  actual  plot  con- 
firms the  reasonableness  of  Equation  (1.1),  including  the  values  of  its 
coefficients. 

In  conclusion,  for  the  situation  considered  here  the  transfer  func- 
tion model  is  more  effective  as  a tool  that  gives  us  some  insight  into  the 
manner  in  which  the  system  operates,  rather  than  as  a tool  that  can  give 
us  reliable  forecasts  of  future  running  times. 

The  remainder  of  this  paper  is  devoted  to  a discussion  of  the  perti- 
nent details  that  lead  us  to  our  model.  In  Section  2,  by  way  of  presenting 
some  aspects  of  transfer  function  model  building,  we  also  introduce  some 
terminology  and  notation.  In  Section  3,  we  present  an  analysis  of  our  data. 

In  what  follows,  we  require  the  reader  to  have  some  familiarity 
with  the  material  In  Box  and  Jenkins  (1976)  and  with  that  in  Jenkins  and 
Watts  (1968). 


T 
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2 . Transfer  Function  Models  and  Their  Estlinatlon 

Univariate  transfer  function  models  as  described  by  Box  and  Jenkins 
(197b)  are  models  that  specify  the  stochastic  interrelationships  between 
two  time  series.  They  are  more  general  than  regression  models  with  lag 
structures  on  predetermined  variables,  in  that  the  dependent  variable  can 
also  have  a lag  structure.  In  addition,  the  transfer  function  models  can 
have  a superimposed  error  structure  which  may  be  of  a very  general  nature. 
Engineers  often  refer  to  error  with  the  term  "noise,"  and  "white  noise" 
refers  to  errors  that  are  independent  and  Identically  distributed. 

There  are  two  equivalent  representations  of  a univariate  transfer 
function  model.  One  is  the  infinite  or  reduced  form,  and  the  other  is  the 
finite  form.  In  the  infinite  form,  the  output  series,  say  (In  our 

case  the  running  time),  is  explicitly  represented  as  a function  of  the  input 
series  (in  our  case  the  down  time)  and  its  lagged  (previous)  values; 

that  is. 


"^t  ~ ''l^t-1  ''2^t-2  •••  ^t  ’ 


(2.1) 


where  the  constants  Vq,v^,...  are  called  the  impulse  response  weights. 

In  cases  where  there  is  no  Immediate  response,  one  or  more  of  the  initial 
v's  , say  '^o’'^l’ ' ' " ’'^b-1  ’ equal  to  zero.  The  process  represents 

noise,  which  is  assumed  to  be  independent  of  the  level  of  the  input  series, 
but  is  additive  with  respect  to  the  influence  of  the  input;  can  have 

any  general  structure. 

It  can  be  shown  that  an  equivalent  representation  of  the  model  given 
by  Equation  (2.1)  is  the  following  finite  form: 

t 1 t-1  r t-r  0 t-b  1 t-b-1  s t-b-s  t 

(2.2) 

where  the  6's  , the  o)' s , and  b are  unknown  constants.  The  constant  b 
associated  with  the  leading  indicator  series  X^  Indicates  which  of  the 

previous  values  of  X^  affect  the  present  Y^  . In  our  application,  the 
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value  of  b represents  the  number  of  previous  maintenance  times  affect- 
ing the  present  running  time. 

A first  step  towards  estimating  the  transfer  function  model  is  a 
tentative  identification  of  the  values  of  r , b , and  s . This  can  be 
accomplished  by  an  examination  of  the  estimated  impulse  response  weights 
Vj^  , k=0,l,2,...  . A plot  of  Vj^  versus  k is  known  as  the  impulse 

response  function. 

There  are  two  general  approaches  for  obtaining  the  Impulse  response 
function.  The  first  one,  outlined  by  Box  and  Jenkins  (1976,  p.  379),  is 
based  on  a "prewhitening"  of  the  input  series.  Prewhitening  the  input 
series  means  fitting  a time  series  model  to  the  series  such  that 

the  residuals  from  the  model,  say  , are  independent  and  Identically 

distributed  random  variables  with  mean  zero  and  a constant  variance. 

When  the  v^^  are  estimated  using  the  prewhitening  of  the  input  series 

procedure,  their  neighboring  values  tend  to  be  correlated.  Thus  the 
graph  of  the  impulse  response  function  tends  to  be  misleading.  This 
ultimately  affects  our  ability  to  obtain  a realistic  transfer  function 
model.  We  are  therefore  Interested  in  considering  an  alternate  approach 
for  estimating  the  impulse  response  weights. 

The  second  approach  for  estimating  the  impulse  response  weights 
involves  the  use  of  "cross  spectral  analysis."  Such  an  approach  removes 
the  difficulties  associated  with  the  problem  of  the  correlated  estimates 
of  V|^  , and  also  provides  us  with  some  additional  insight  into  the  na- 
ture of  the  dependencies  between  the  input  and  the  output  series.  These 
are  illustrated  at  the  end  of  Section  3. 


Once  the  impulse  response  function  is  obtained,  we  can  isolate  the 
noise  series,  , by  using  Equation  (2.1).  Specifically,  we  estimate 

the  noise  series  by 


^t  ~ ''l^-l  " ’^2^t-2 
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A knowledge  of  , plus  a knowledge  of  the  tentative  values  of 

r , b , and  s , helps  us  to  estimate  the  parameters  of  the  transfer  func- 
tion model.  Equation  (2.2).  One  way  of  accomplishing  this  is  by  using  the 
TIMES  program  package  described  by  Willie  (1977). 

The  adequacy  of  the  proposed  model  can  be  checked  by  an  analysis 
of  the  residues  from  the  model.  The  details  of  such  an  analysis  are 
given  in  Box  and  Jenkins  (1976,  p.  392). 

3.  Analysis  of  Running  Time  and  Down  Time  Data 

Data  on  the  operating  history  of  nuclear  reactors  are  generally 
available  showing  dates  on  which  the  reactors  ceased  operation  and  the 
duration  of  the  stoppage.  Among  other  facts,  the  reasons  for  the  stop- 
page are  also  given.  Stoppages  are  categorized  according  to  whether  they 
were  scheduled  or  forced.  In  addition  to  this,  there  is  a further  break- 
down indicating  whether  the  stoppage  was  due  to  equipment  failure,  testing, 
refueling,  regulatory  reasons,  operator  training,  administrative  reasons, 
operational  error,  or  other  causes. 

3.1  Screening  the  data 

The  dat.a  that  were  given  to  us  described  the  Robinson  Power  Plant's 
operating  history  from  June  1974  through  April  1976.  These  data  did  not 
contain  any  stoppages  due  to  administrative  reasons,  operational  error,  or 
other  casues;  they  contained  one  stoppage  for  regulatory  reasons  and  one 
stoppage  for  refueling.  In  one  instance,  the  data  contained  an  unrealistic 
combination  for  the  cause  of  stoppage — a scheduled  failure.  In  this  case 
we  used  our  discretion  to  alter  it  to  a forced  failure.  Whenever  there 
were  stoppages  due  to  operator  training,  these  ware  treated  as  running  times 
rather  than  as  down  times.  This  was  done  for  two  reasons.  First,  the  du- 
ration of  each  stoppage  was  very  short  (on  the  average  about  half  an  hour); 
second,  we  would  like  to  concentrate  on  those  down  times  that  pertain  to 
the  physical  operation  of  the  system  rather  than  on  those  external  to  it. 
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The  single  stoppage  due  to  a forced  regulatory  restriction  was 
for  a period  of  3.67  hours,  and  since  it  immediately  followed  a forced 
equipment  failure  of  15.52  hours,  it  was  combined  with  the  equipment  fail- 
ure stoppage. 

Refueling  the  reactor  takes  place  annually  and  is  generally  of  a 
very  long  duration.  In  our  data,  we  had  only  one  stoppage  for  scheduled 
refueling,  and  it  was  of  960  hours'  duration.  Since  the  duration  of  this 
stoppage  is  out  of  line  with  the  duration  of  the  other  stoppages  (see 
Table  1.1),  it  was  excluded  from  consideration. 

We  remark  here  that  any  time  a stoppage  (running  time),  say  ^^t^  ’ 

was  excluded  from  consideration,  its  corresponding  running  time  (stoppage) 

^t  ^^t^  also  excluded.  This  is  to  ensure  that  no  bias  is  introduced 

into  the  relationship  between  the  two  variables  of  interest  because  of  the 
elimination  of  observations  of  either  one. 

Another  convention  followed  in  our  analysis  arises  from  the  fact 
that  the  raw  data  show  the  dates  and  the  duration  (in  hours)  of  each  stop- 
page, rather  than  the  actual  time  of  stoppage.  We  assume  that  each  down 
time  commenced  at  0000  hours  (unless  in  some  rare  instances  there  is  a 
second  stoppage  occurring  during  the  same  day).  Whenever  two  or  more 
breakdowns  occurred  during  the  same  day,  they  were  combined  into  one  down 
time  period  and  the  intervening  operating  period  was  ignored. 

Because  of  the  paucity  of  data,  we  chose  not  to  distinguish  be- 
tween stoppages  due  to  equipment  failure  and  those  due  to  testing.  This 
is  reasonable  because  whenever  there  is  a forced  equipment  failure,  main- 
tenance and  test  actions  on  other  (nonfailed)  components  are  routine. 

Thus,  in  practice  it  is  difficult  to  differentiate  clearly  between  the 
consequences  of  equipment  failure  and  those  of  testing.  This  strategy 
was  suggested  by  some  staff  members  of  the  Probabilistic  Analysis  Staff 
at  NRC. 


The  preliminary  screening  and  examination  described  above  gives 
us  a series  of  values  of  the  down  times  (in  days)  and  a series  of 
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values  of  the  corresponding  runi'ing  times  (In  days).  However,  as 

we  shall  soon  see,  some  further  screening  is  necessary. 

The  next  step  in  our  analysis  involved  prewhitening  the  series. 

This  turned  out  to  be  quite  a frustrating  endeavor,  since  no  simple  uni- 
variate time  series  model  of  the  Box- Jenkins  type  seemed  to  provide  a 
reasonable  fit.  The  difficulty  turned  out  to  have  been  caused  by  two 
unusually  large  down  times  due  to  scheduled  testing  and  forced  failure 
of  18.08  and  25.43  days,  respectively.  These  were  incompatible  with  the 
other  down  times  (see  Table  1.1),  and  thus  defied  the  use  of  a simple  model 
as  a prewhltenlng  transformation.  Perhaps  a time  series  model  with  an 
indicator  variable  [such  as  those  used  in  the  "Intervention  analysis"  of 
Box  and  Tiao  (1975)]  might  have  been  adequate  for  these  and  for  the  re- 
fueling stoppage,  but  this  was  not  attempted.  In  the  Interest  of  expe- 
diency, it  was  preferable  to  eliminate  the  two  large  X^'s  and  their 

corresponding  Y^’s  . Thus,  in  effect,  some  data  screening  was  done  during 

the  prewhitening  phase.  Table  1.1  presents  the  28  screened  values  of  the 
down  times  and  the  corresponding  running  times  . We  remind  the 

reader  that  the  subscript  t is  a sequential  index  rather  than  an  index 
representing  time.  That  is,  X^  and  Y^  are  not  simultaneously  observed 

in  time;  X^  precedes  Y^  . 

3.2  Transfer  function  model  development 

In  Figu’" 'S  3.1  and  3.2  we  show  plots  of  the  autocorrelation  func- 
tions of  the  (screened)  down  times  X^  and  the  (screened)  running  times 

Yj_  , respectively.  Based  upon  these  plots  we  are  able  to  conclude  that 

the  two  time  series  can  be  treated  as  stationary  (see  Box  and  Jenkins 
(1976),  p.  174). 

In  Figure  3.3  we  show  a plot  of  the  logarithm  of  the  smoothed  es- 
timator of  the  power  spectrum  of  the  screened  down  times  X^  , as  well  as 
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Figure  3.3  smoothed  estimator  of  the  spectrum  of  down  times  versus 

frequency,  using  rectangular  window  of  band  width  .098. 
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the  approximate  95%  confidence  limits.  The  smoothing  was  performed 
using  a rectangular  window  of  band  width  .098.  As  a matter  of  fact, 
all  the  smoothing  that  is  discussed  in  this  paper  was  performed  using 
a rectangular  window  of  band  width  .098.  The  power  spectrum  curve  shows 
us  how  the  variance  of  the  time  series  is  distributed  with  frequency. 

For  a detailed  understanding  of  the  power  spectrum,  its  smoothing,  and 
the  band  width  of  a smoothing  window,  we  refer  the  reader  to  Chapter  6 
of  Jenkins  and  Watts  (1968). 

Our  next  step  involves  the  determination  of  a suitable  prewhitening 
transformation  for  the  down  times  . Based  upon  the  several  models  that 

we  attempted,  we  conclude  that  a moving  average  process  of  order  3 best 
describes  the  X^  series.  Specifically,  we  find  that 

X^  - 0.7963  = + 0.117a^_^  - 0.189a^_2  - 0.133a^_2  , (3.1) 

where  .7963  is  the  mean  of  the  X^  series.  The  represent  the  re- 

siduals when  a moving  average  process  of  order  3 is  fitted  to  the  X^ 

series.  If  the  prewhitening  transformation  given  by  Equation  (3.1)  is 
correct,  then  the  a^.  will  be  Independently  and  identically  distributed 

with  a constant  mean  and  variance. 

In  Table  1.1  we  give  the  values  of  the  a^'s  . In  order  to  verify 

the  appropriateness  of  the  model  given  by  Equation  (3.1),  we  plot  the  esti- 
mated autocorrelation  function  and  the  estimated  power  spectrum  of  the 

series.  These  plots  are  given  in  Figures  3.4  and  3.5,  respectively.  We 
remark  that  the  plot  of  the  estimated  power  spectrum  of  the  series 

given  in  Figure  3.5  is  relatively  constant  as  compared  to  the  plot  of 
the  estimated  power  spectrum  of  the  X^  series  given  in  Figure  3.3. 

This  is  because  the  effect  of  prewhltenlng  is  to  remove  the  dependencies 
among  the  ^nd  give  us  a set  of  Independent  . Figures  3.4  and 

3.5  confirm  the  appropriateness  of  the  prewhltenlng  transformation  given 
by  Equation  (3.1). 
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Following  Box  and  Jenkins  (1976,  p.  380),  we  next  apply  the  same 
prewhitening  transformation  [i.e.,  the  one  given  by  Equation  (3.1)]  to 
the  running  times  and  obtain  the  as  residuals.  In  Table  1.1 

we  give  the  values  of  3^  under  the  heading  "Prewhitened  Running  Times." 

In  Figure  3.6  we  show  a plot  of  the  smoothed  estimator  of  the  power  spec- 
trum of  . We  remark  that  except  at  the  very  low  frequencies,  the 

plot  of  the  estimated  power  spectrum  of  the  fairly  constant. 

Thus  it  appears  that  the  prewhitening  transformation  given  by  Equation 
(3.1),  when  applied  to  the  running  times  , also  yields  a sequence  of 

independent  and  identically  distributed  random  variables  3^.  • 

Our  next  step  is  to  obtain  the  cross  correlation  between  the 

and  P at  lags  k , k=0,l,2,...  . If  s and  s„  denote  the  estimated 

t Ot  p 

standard  deviations  of  the  and  the  3^  series,  respectively,  and  if 

r^g(k)  denotes  the  estimated  cross  correlation  between  the  and  the 

3j.  at  lag  k , then  v^^  , an  estimate  of  Vj^  , is 

\ " 'T  ' ...  , 

a 

(see  Box  and  Jenkins  (1976),  p.  380). 

In  Figure  3.7  we  show  a plot  of  the  impulse  response  function; 
that  Is,  a plot  of  v^^  versus  k , k=0,l,2,...  . We  remark  that  in 

Figure  3.7,  the  value  v^  is  significantly  larger  than  the  other  values 

of  V|^  , k=l,2, ...  . 

Because  the  neighboring  values  of  the  Vj^  in  Figure  3.7  tend  to 

be  correlated,  we  also  obtained  the  Impulse  response  function  using  the 
"cross  spectrum"  between  the  and  the  Y^  series  (see  Jenkins  and 

Watts  (1968),  p.  A24).  In  Figure  3.8  we  show  a plot  of  the  Impulse 
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,1  m.1  jor  ]ioiiU  of  our  i-om  Iiisicins. 


» 'c  '(;  ' i;,'  r * . w 

*.  This  is,  of  course. 


We  now  estim.ite  the  noist'  series  usinp  Kquat  ion  (7.3';.  An 

jiispeclion  of  the  e,stim.)led  autocorrelation  function  of  the  estimated 
noise  series  , .ind  ,i  "portmani  eau  l.ick  of  fit  test"  (l.hP  with  8 

dep.rees  of  freedom)  (see  liox  and  denkins  (1976),  p.  290),  lead  us  to 
conclude  that  the  .ire  independent  and  identical Iv  distributed.  In 

.'iddition  to  the  above,  we  sliow  in  Figure  3.9  a plot  of  the  smoothed  es- 
I imator  of  the  power  spectrum  of  the  . Here  again,  except  at  the  verv 

low  frequencies,  the  estimated  power  spectrum  of  the  series  is  fairly 

const. mt.  Thus  tlie  estimated  noise  series  N can  be  described  rcason- 

t 

.ably  Well  by  .a  white  noise  [iroeess. 

A final  step  in  the  analysis  involves  the  fitting  of  a transfer 
function  model  to  the  running  times  . This  w.as  accomplished  by  using 

the  I'lMFS  program  package.  Of  the  several  models  that  were  attempted, 
the  model 

Yj  - 10.630  = 6.6X^  - U.5'3.X^_^  (3.2) 

appears  to  be  the  best;  10.630  is  the  mean  of  the  Y^  series. 


In  order  to  verify  the  reasonableness  of  the  model,  two  diagnostic 
checks  were  suggested  in  Box  and  Jenkins  (1976).  One  depends  on  the  auto- 

A 

correl.it  ion  of  the  .ind  the  other  depends  on  the  estimated  cross  cor- 

A 

relation  between  the  and  the  . For  both  cases  a portmanteau  lack 

of  fit  test  was  used.  In  the  former  case,  the  test  statistic  is  1.69  with 
8 degrees  of  freedom,  whereas  in  the  latter  case  the  test  statistic  is 
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.69H  wiLti  6 df>>r(.'cs  ol  I rot’dom.  These  values  support  the  reasonableness 
ol  1 In-  proposed  transfer  function  model. 


As  stated  enrlii'r,  a cross  spectral  analysis  of  the  and  tiie 

series  can  f>ive  us  further  insiglit  into  the  nature  of  the  dependen- 

l ies  hi-tween  thi‘  two  series.  For  example,  G(f)  , the  "sain"  at  frequen- 
i-y  f , hi-haves  like  tiu'  regression  coefficient  in  a linear  regression 
model  hi'tween  the  output  and  the  input  at  frequency  f (see  Jenkins  and 
Watts  (1968),  p.  332).  In  Figure  3.10  we  plot  the  gain  of  the  running 
limes  on  tlie  down  times  at  various  frequencies.  Another  important  func- 
l ion  in  cross  spectral  analysis  is  the  "squared  coherence"  between  the 
input  and  the  output  at  frequency  f . This  quantity  measures  the  cor- 
relation between  the  sinusoidal  component  of  Y and  that  of  X at 

t t 

frequency  f . The  square  coherence  is  also  in  some  sense  a measure  of 
the  proportion  of  information  in  the  series  that  is  attributable 


to  the  X series, 
t 


For  more  information  on  the  coherence  and  the  co- 


herence spectrum,  we  refer  the  reader  to  Jenkins  and  Watts  (1968,  p. 

332).  In  Figure  3.11  we  show  a plot  of  the  coherence  for  the  running  time 
and  down  time  data.  Figures  3.10  and  3.11  give  us  some  additional  assur- 
ance on  the  dependence  of  the  running  times  on  the  down  times. 


4 . Summary  and  Conclusions 

In  the  foregoing  analysis  we  have  demonstrated  the  use  of  time 
series  analysis  methodology  for  studying  the  interrelationships  between 
the  maintenance  times  and  the  running  times  of  a complex  system.  Our 
analysis  enables  a decision  maker  to  assess  the  impact  of  his  maintenance 
policies  on  running  times,  or  to  influence  the  operating  times  by  managing 
the  maintenance  times.  In  addition,  given  a down  time  a decision  maker 
can,  to  a tirnitrd  extemi  , forecast  the  next  running  time.  This  type  of 
information  may  be  very  valuable,  especially  for  large  and  complex 
syst(?ms . 
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Our  analysis  can  be  criticized  on  the  grounds  that  it  is  based  on 
an  i nsiiT r ic ient  amount  oi  data.  We  hope  that  this  criticism  can  be  over- 
looked in  the  liglit  of  the  fact  that  our  approacli  is  to  he  viewed  as  a 
prototype  f('r  tlu'  analysis  of  reliabilitv  data  involving  two  interrelated 
sources  of  data. 
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